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CN ■ Abstract 

C^) | Acoustic emission before the failure of heterogeneous materials is studied as a func- 

tion of applied hydrostatic pressure. A formula for the energy release is suggested, 
which is valid in the whole diapason of pressures, from zero to the critical pressure 
of rupture. This formula is obtained by employing the extrapolation technique of 
the self-similar approximation theory. The result is fitted to experiment in order to 
demonstrate the correct general behaviour of the obtained expression for the energy 
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1 Introduction 

There exists a widespread understanding in scientific community that the moment of rupture 
in heterogeneous materials is somewhat similar to a critical point. The fracturing process is 
then a kind of a second-order phase transition [2,4]. This underlies the intuitive assumption 
which compares global failure in disordered materials at mesoscopic scale with percolation 
at microscopic scale. So, as percolation, global failure is assumed to be a critical phe- 
nomenon. The fracturing process is then classified as a self-organised random irreversible 
process. Extensive numerical simulations and experimental measurements have confirmed 
this conclusion. This concerns the global failure of materials as well as earthquakes [9], since 
an earthquake can also be considered as a global failure of a large material mass. The classi- 
cal power law, defining the energy release rate, has been verified in the critical region, close 
to the time of either a large earthquake or a global materials failure. The related critical 
exponent is, in general, a complex number, which is due to taking into account long-range 
elastic interactions transported by stress fields around defects and cracks. The imaginary 
part of the exponent gives rise to the so-called log-periodic corrections which have been 
identified quite early in renormalization group solutions for critical phase transitions [7,10]. 
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Experimental results of acoustic emission measurements and seismograms unambiguously 
revealed the existence of oscillatory corrections before materials fractures and earthquakes, 
respectively [1,10]. If one would know the general law describing this acoustic emission, one 
would be able to make an early prediction of the global failure in heterogeneous materials. 
It is a principal, not yet solved problem, to find such a general law, which would be valid not 
solely in the asymptotic vicinity of the critical point, but in the whole interval of pressures 
or times, both near the critical point of failure as well as far from it. And this is the main 
aim of the present paper to suggest a mathematically grounded derivation for such a law 
and to confront it to available experimental data [5] . For this purpose, we employ the recent 
progress in the self-similar approximation theory [12-18] which has been successfully applied 
to the theory of critical phenomena and to time series forecasting [19-21]. We would like to 
stress here that this theory is general and can be used for any given series of experimental 
data, whether this concerns phase transitions, financial markets, or any other series. 

2 Modelling 

Let us be interested in the behaviour of the energy release E(jp) as a function of pressure p, 
considered between p !=s 0, i.e. at the initial stage when the very early damages occur, and up 
to the point p = p c of global failure. Note that p refers to fracturing under applied spatially 
uniform loading, as in experiments where water pressure inside a spherical tank increases 
up to the material failure. Experimentally, the applied pressure is usually increased linearly 
with time. The behaviour of E(p) at the vicinity of the critical point of rupture is 

E{p) w E{p c ) + A{p c - p) a + B{p c - p) a cos[cu ln(p c - p) + <p] , (1) 

as p — > pc, where A, B, a, u and ip are parameters (a and u are the real and imaginary parts, 
respectively, of the so-called critical complex exponent). The problem is how to find the 
behaviour of E(p) in the whole range < p < p c l Here we consider the case of one variable, 
pressure. A more general loading, involving not one variable, like p, but several variables 
can also be treated by the theory. However, before being involved in such generalisations, 
we would like to show on a simpler example how to accomplish the extrapolation. Our 
main aim in this letter is to demonstrate the general possibility of realising such an accurate 
extrapolation for the energy release. 

Let us transform Eq. (1) by introducing the dimensionless variable x = (p c — p)/p c , 
(0 < x < 1) and the reduced energy release f(x) = E{p)/E{p c ) which is a dimensionless 
function. Then, expansion (1) takes the form 

f(x) ps 1 + a(x) x a , a(x) = A[l + cos(c<jln x + (3)} (2) 

whose parameters /3, A, and fi can be easily expressed through p c , E(p c ), A, B, uj, and (p. 
Note that function a(x) cannot be expanded in powers of x near zero. Therefore, expansion 
(2) can be treated as a generalised asymptotic series, as defined by Poincare [8], with a(x) 
considered as a coefficient. Employing the self-similar approximation theory [12-21], the 
asymptotic series (2) can be extrapolated to 

f*(x) = exp[c(x) x a ] , c(x) = [1 + yucos(o;lna; + ft)] r , (3) 
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where r has to be defined from an optimization or boundary condition. As is clear, the 
boundary condition E(0) ~ is not convenient here, leading to large errors in defining r 
because of the physical impossibility to determine the very early precursors with a high 
accuracy. A more judicious choice is to consider a global physical quantity. In the same way 
as treating E{p) as the integrated energy release rate, we may define F{p) as the integral 
energy release Jq E(q)dq. Given p , the parameter r then can be defined from the sum 
optimization rule 

L~™*=j5& (»-°. »-*.-<»■ w 

In general, the lower limit in the integral (4) can pertain to the whole interval [0,1]. In reality, 
the trustful experimental data are available only below p c . One has to choose the data at the 
highest available pressure, so that the lower limit in the sum rule (4) be sufficiently small. 
In this way, the self-similar approximant (3) becomes 

f*(x) =exp{[l + /icos(cc>lna; + / 3)] rx a } (5) 

with r given by the sum rule (4). Note that this formula contains the same number of 
parameters as the initial formula (1). Returning, with the help of relations (1) and (2), to 
the energy release, we come to the formula 

E*{p) = E(p c ) f*(x) , (6) 

which extrapolates E(p) for the whole interval < p < p c . It is worth stressing that the law 
(6) is obtained as a self-similar extrapolation of the acoustic emission signals observed at 
the vicinity of the rupture. Another possibility would be to base this kind of extrapolation, 
starting from the signals existing at the initial stage of the process, being yet very far from the 
rupture. This latter way was used by Gluzman et al. [3], who started with a given polynomial 
expansion at the early time of acoustic emissions, treating the power law behaviour close to 
rupture as a boundary condition. To our mind, this opposite approach has two weak points: 
First of all, the early acoustical precursors, as is known, are very small, being embodied in 
the acoustical noise, because of which they can be hardly measured with a good accuracy. 
Second, it looks difficult or even impossible to extract information on log-periodic corrections 
from the early acoustic signals. Another phenomenological expression for the energy release 
in the noncritical region has been suggested by Sornette and Andersen [10] whose arguments 
were based on the existence of a scaling of the macroscopic elastic modulus and the elastic 
energy release rate in the thermal fuse model [6]. The following phenomenological form has 
been proposed [5] 

E(p) « E(p c )A (tanh ^-^) ° + B (tanh ^) ° cos jw (tanh ^-^) + <p] , (7) 

which is a pure power law, like Eq. (1), in the critical region close to rupture and exponen- 
tially relaxes in the noncritical region far from rupture, where only a few damages occur. 
However, the suggested form (7) possesses an unphysical behaviour, being negative in a wide 
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range of its variable. At the same time, our self-similar formula (5) yields the energy release 

(6) that is always positive. In Fig. 1, we compare the behaviour of E{p) , in the dimen- 
sionless form, for both Eqs. (6) and (7), with the parameters chosen so that to have similar 
variations when approaching the critical point. In this figure, one can clearly see that Eq. 

(7) possesses unphysical negative values for the major part of the interval [0,1], while Eq. 
(6) is everywhere positive. 

3 Comparison with experimental data 

In order to fit our self-similar formula to experimental results of acoustic emission measure- 
ments, we refer to the most accurate data, available nowadays, that is those of Anifrani's team 
[1] at Aerospatiale-Matra Inc., which are reported by Johansen and Sornette [5]. In these 
measurements, the acoustic emission has been recorded during the pressurisation of spher- 
ical tanks of kevlar or carbon fibres pre-impregnated in a resin matrix wrapped up around 
a thin metallic liner (steel or titanium) fabricated and instrumented by Aerospatiale-Matra 
Inc. It has been found that the seven acoustic emission recordings of seven pressure tanks, 
which were brought to rupture, exhibit a clear acceleration in agreement with a power law 
divergence as expected by the critical point theory. A strong evidence of oscillatory correc- 
tions, forming the intermittent succession of bursts of acoustic emission, when approaching 
the rupture, has been clearly identified. 

In Fig. 2, we present the energy release rate, which is the derivative of our formula (6) 
for the cumulative energy release. The result is obained by fitting to the experimental curve 
reported at the right bottom in figure 1 of the paper by Johansen and Sornette (2000). These 
experimental data explicitly show an oscillatory behaviour of the energy release rate before 
the rupture. The best fit for the five parameters a, f3, /i, r, and u have been obtained for 
0.7, 1.9, -0.045, 7.5 and 15, respectively, under the given pressure of rupture at 673 bars. 
The parameter r satisfies Eq. (4). The upper pressure po ~ 665 bars, measured before the 
rupture, corresponds to the value yo ~ 0.0119, which is sufficiently small, in agreement with 
condition (4). 

4 Conclusion 

The parameters can change for different materials as well as for different experimental se- 
tups. However, the three key physical parameters: a, u, and r should be universal, not 
essentially depending on the details of the external stress in the critical region. To prove 
this assumption, it is necessary to accomplish a set of experiments for appropriate materials. 
Selecting particular materials, we should keep in mind that: (i) the width of the critical 
region is linearly proportional to the strength of disorder, (ii) long-range elastic interactions 
are responsible for log-periodic corrections, (iii) the more brittle is the fracturing process, 
the stronger are the elastic interactions and stronger the energy release fluctuates, (iv) the 
fluctuations in the energy release are also stronger when a high value of a material charac- 
teristic is alternated with a low value of that characteristic. Thus, elastic-porous materials 
are expected to be the best choice. 
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Figure Captions 



Fig. 1. General form of the dimensionless energy release corresponding to Eq. (6), with 
a = 0.7, uj = 6.9 and r = 6.6 (continuous line) and to Eq. (7), with a = 0.7, u> — 8, and 
t = 0.8 (broken line). 



Fig. 2. Dimensionless energy release rate from Eq. (5), fitted to the experimental data 
reported by Johansen and Sornette (2000). 
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